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Abstract. For axially symmetric flow of dissipationless inhomogeneous fluid onto a non 
rotating astrophysical black hole under the influence of a generalized pseudo-Schwarzschild 
gravitational potential, we investigate the influence of the background flow configuration on 
determining the salient features of the corresponding acoustic geometry. The acoustic horizon 
for the aforementioned flow structure has been located and the corresponding acoustic surface 
gravity k as well as the associated analogue Hawking temperature Tah has been calculated 
analytically. The dependence of k on the flow geometry as well as on the nature of the back 
ground black hole space time (manifested through the nature of the pseudo-Schwarzschild 
potential used) has been discussed. Dependence of the value of k on various initial boundary 
conditions governing the dynamic and the thermodynamic properties of the background fluid 
flow has also been studied. 
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1 Introduction 

Contemporary works in the field of analogue gravity phenomena have attracted significant at- 
tention in the community [1-4]. Proper equivalence has been established between the physics 
of the propagating acoustic (and acoustic type) perturbations embedded in an inhomogeneous 
dynamical fluid system, and some kinematic features of the general theory of relativity. Such 
formalism has opened up the possibility of simulating various important features of the black 
hole space time within the laboratory set up. 

Conventional works in this field, however, concentrates on the physical systems for which 
gravity like effects are realized as emergent phenomena. Such systems do not usually contain 
any source that produces active gravitational field in any form. In recent years, though, 
attempts have been made to study the analogue effects in strong gravity environment [5-13]. 
The uniqueness of such systems lies in the fact that those are the only analogue models 
studied so far that simultaneously contain both kind of horizons, the gravitational as well as 
the acoustic - allowing one to go for a close comparison between the actual and the analogue 
Hawking effect. 

Till date, analogue effects in the axisymmetrically accreting black hole systems has been 
studied for the flow structure assumed to be in hydrostatic equilibrium in the vertical di- 
rection. Two other configurations for the axially symmetric black hole accretion are also 
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possible. Acccretion flow in conical equilibrium [14] and 'flat disk' kind of flow with con- 
stant flow thickness have also been studied in the literature (see, e.g., [15] for detail review). 
For pseudo-Schwarzschild gravitational potential x , the aforementioned three different flow 
structures have quite recently been studied in great detail ([18]). to reveal various critical 
phenomena and stability properties of such flows. 

In our present work, we study the analogue effects in accreting black hole systems under 
the influence of pseudo-Schwarzschild potentials for three different flow configurations - flow 
in conical and in vertical equilibrium and for constant thickness (height) flow. The main 
motivation behind this work is to study the dependence of the salient features of the acoustic 
geometry on the background geometry realized through the aforementioned three different 
flow configurations. To accomplish such task, we will calculate the acoustic surface gravity 
for same set of accretion parameters but for all possible flow configurations in a pseudo- 
Schwarzschild potentials. 

2 Acoustic surface gravity for classical analogue systems 

Classical analogue gravity systems (alternatively, the classical 'black hole analogues') are 
fluid dynamical analogue of black holes in general relativity. Such analogue may occur when 
a small linear perturbation propagates through a dissipation less inhomogeneous barotropic 
transonic fluid at finite temperature. The corresponding acoustic metric, which specifies 
the geometry in which the perturbation propagates, may be constructed in terms of the 
flow variables defining the unperturbed background continuum. The transonic surface acts 
as acoustic horizon - a null hypersurface with acoustic null geodesies, the phonons, as its 
generators. The acoustic black hole horizon, which resembles the black hole event horizons 
in many ways, may form at the regular transonic point of the fluid, whereas an acoustic white 
hole horizon may be formed at the hypersurface where the fluid makes a discontinuous sonic 
transition, e.g., through a stationary shock [8, 19]. 

In his pioneering work Unruh [20] introduced the concept of acoustic geometry inside 
a supersonic fluid and demonstrated that an analogue surface gravity n may be associated 
with an acoustic black hole type event horizon, and one of the most interesting aspects 
of the acoustic horizon is to emit the Hawking type radiation of thermal phonons. Such 
acoustic Hawking radiation may be characterized by a analogue Hawking 2 temperature Tah = 

In Unruh's original approach, the acoustic surface gravity could be associated with the 
component of the bulk velocity of the flow normal to the acoustic horizon u± and the speed 
of propagation of the acoustic perturbation c s as 



where J-J- = rfd^ represents the space derivative taken along the normal to the acoustic 
horizon, and every quantities in the eq. (2.1) has been evaluated at the location of the 
acoustic horizon r^. 

order to optimize between the easy handling of the Newtonian framework of gravity and more rigorous 
and non-tractable complete general relativist ic description of the strong gravity space time, four different 
'modified' Newtonian 'black hole potentials' have been introduced in the literature which are commonly 
known as pseudo-Schwarzschild potentials, see, e.g., [16] , and [17] for further detail about such potentials. 
2 Hereafter the phrases 'acoustic' and 'analogue' will be used synonymously for the sake of brevity. 
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In eq. (2.1), however, the sound speed was assumed to be a position independent constant. 
Unruh's work was followed by several other important contributions ([21-24] , to name a 
few) 3 Visser [23] implemented the additional contribution due to the position dependent 
sound speed and obtained a modified expression for the surface gravity 
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The acoustic horizon is a surface defined by the equation [23] 
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(2.2) 
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for the stationary background flow configuration. Equation (2.3) basically states that the 
acoustic horizons is a transonic surface. The supersonic region of a transonic flow defines the 
acoustic ergo region. 

The concept of acoustic geometry has been extended to a relativistic fluid flow in a 
general background space time [26]. For an ideal barotropic fluid, the relativistic Euler 
equation and the equation of continuity obtained from the energy momentum conservation 
can be linearized in order to obtain the wave equation for the propagating perturbation in 
analogue curved space time with the corresponding acoustic metric. The generalized form of 
the acoustic surface gravity turns out to be 



(1 



drj 



(u± - c s 
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rh 



where x M is the Killing field which is null on the corresponding acoustic horizon. The alge- 
braic expression corresponding to the v / 3v 7r \> ma y thus be evaluated once the background 
stationary metric governing the fluid flow as well as the propagation of the perturbation in a 
specified geometry with well posed boundary conditions are realized. It is worth mentioning 
that the generalized form for k as defined in eq. (2.4) can further be reduced to its New- 
tonian/semi Newtonian counterpart depending on the nature of the gravitational potential 
describing the background fluid motion. 



3 Acoustic surface gravity for axisymmetric black hole accretion 

From (2.2) and (2.4) it is clear that, in order to find the acoustic surface gravity n for the 
Newtonian as well as for the general relativistic acoustic geometry, it is sufficient to calculate 
the location of the acoustic horizon, the sound speed c s of the small linear perturbation 
and its normal space gradient dc s /drj, as well as the the normal (to the acoustic horizon) 
component of the flow velocity u± and its normal space gradient du±/drj, evaluated on the 
acoustic horizon. 

For transonic accretion onto astrophysical black hole, one thus needs to consider the 
Euler equation and the equation of continuity for a specific symmetry of the problem. The 
Euler and the continuity equation may then be linearized (for a general linearization scheme 
and its application to the study of axisymmetric black hole accretion for three different flow 

3 It is interesting to note that the concept of the acoustic geometry inside a transonic fluid was first realized 
by [25] while studying the stability properties of the relativistic spherical accretion onto astrophysical black 
holes. 



-3- 



geometries see [18]) in order to construct the corresponding acoustic metric to specify the 
relevant acoustic geometry. The structure of the stationary background fluid flow may be 
provided by the stationary solution of the Euler and the continuity equations. For a certain 
set of values of the initial conditions describing the transonic accretion flow governed by 
certain barotropic equation of state, it may be possible to find the location of the saddle type 
transonic point, which is identical to the acoustic black hole horizon (see [7, 8, 10, 26] ). The 
quantities [u±, c s , du±/dr, ^ c s/^ r ] rtransonic subject to a gravitational field which governs the 
accretion, may then be evaluated to estimate the acoustic surface gravity using eq. (2.4). 
The quantity y/x^Xn evaluated at the acoustic horizon is a function of = r tra nsonic (as will 
be demonstrated in subsequent paragraphs) depending on initial conditions. 

In general, we consider the equatorial slice of the axisymmetric gravitating accretion of 
hydrodynamic fluid onto non rotating black holes. The gravitational field of the black hole 
is assumed to be described by certain pseudo-Schwarzschild potentials. The Schwarzschild 
radius r g = 2GMbr/c 2 is used to scale the radial distance, whereas all velocities involved 
are scaled by the velocity of light in vacuum c, G = c = 1 has been used. Accretion is 
assumed to possess finite radial velocity u commonly known as the 'advective velocity' in 
the accretion literature. Considering v to be the magnitude of the three velocity, u is the 
component of three velocity perpendicular to the set of timelike hypersurfaces {^ v } defined 
by v 2 = constant. The advective velocity u is thus perpendicular to the acoustic horizon 
and hence u is identical with u±. Hereafter we drop the subscript _L in u± and simply use u 
instead. 

The low angular momentum sub-Keplerian advective inviscid flow will be considered 
where the specific flow angular momentum A will be assumed to be a position independent 
constant. Viscous transport of angular momentum will not be taken into account since close 
to the black hole, the infall time scale for the highly supersonic flow is rather small compared 
to the corresponding viscous time scale (see, e.g., [17, 27, 28] and references therein for 
further detail). Also for advective accretion, large radial velocity at a larger distances are the 
consequence of the small rotational energy of the flow [29-31] . The corresponding angular 
velocity £1 may be defined as ([28] and references therein) 

where gij are the metric components. 

The corresponding surface gravity k can now be expressed as 



— g rr 1 — c 2 s \_dr 



{u - c s ) 



(3.2) 

rh 



where X M — C M + ^ M an d the Killing vectors ^ and ^ are the generators of the temporal 
and the axial symmetry group. The norm of the Killing vector Xu may be computed as 

i 

= (g tt + 2ttg t4> + n 2 g4>4> ) 3 = — (3.3) 



where 

^2 2 



£ = 9t<t> - 9ttg<p<p (3.4) 
A 2 = (g„ + 2\g t4> + \ 2 - 



-4- 



In the Newtonian limit 
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where $ is the pseudo potential. Hence the acoustic gravity for axisymmetric black hole 
accretion under the influence of the pseudo-Schwarzscild potential becomes 
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Our main task now boils down to the evaluation of c s , dc s /dr and du/dr on the acoustic 
horizon for the pseudo-Schwarzschild axisymmetric black hole accretion in three different 
flow configurations as mentioned in previous sections. 



4 Generalized transonic accretion in three different flow geometries 

The governing equations describing the dynamics of axially symmetric pseudo-Schwarzschild 
inviscid hydrodynamics accretion are the equation for the conservation of linear momentum 
(the Euler equation): 



9 , v , x d , N 1 dp(r.t) A 2 



where u, p and ^, being the dynamical flow velocity, the fluid density and the pressure, 
respectively, are functions of both r and t. $ may be taken as any one of the following four 
pseudo-Schwarzschild potentials: 

= ~w=Ty (42a> 

1 

2r 

$3 = -1 + (1 - , (4.2c) 

* 4 = -ln(l--J (4.2d) 



$ 2 = _J- 
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The potential $i and $2 have been introduced in [32] and [33] , respectively, whereas $3 
and $4 have been introduced in [16] . 

The mass conservation equation (the continuity equation) can be written as: 

g- t p(r, t) + — [p(r, t)u(r, t)rH] = (4.3) 

The quantity H is the flow thickness which is different in three different flow configurations. 
For the simplest possible flow configuration, the flow thickness (i.e., the disc height) is con- 
stant, and hence H is not a function of the radial distance. In its next variant, the flow can 
have a conical structure (see [14] ) where H is directly proportional to the radial distance 
as H = Ar, where the geometric constant A depends on the solid angle subtended by the 
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flow. For hydrostatic equilibrium in vertical direction [34, 35], the flow thickness can have 
a rather complex dependence on the radial distance and on the local adiabatic sound speed 

[17, 27] defined as c s = (j^J • The barotropic equation of state will in general be used 

to describe the accretion flow in this work. A polytropic equation of state p = Kp 1 will be 
used, whereas the isothermal flow will be governed by the equation p = The quantities 

K,7, kb,T, [i and m# are the entropy per particle, the Boltzmann constant, the isothermal 
flow temperature, the reduced mass and the mass of the Hydrogen atom, respectively. With 
the help of the equation of states, and the specified radial dependence of the flow thickness, 
we can find stationary solutions of the Euler and the continuity equations and draw the Mach 
number versus the radial distance phase portrait for the integral flow solutions to obtain the 
detail information about the location of the acoustic horizon as well as the horizon related 
quantities. 

4.1 Polytropic Accretion 

For polytropic equation of state, the integral solution of the stationary part of the Euler 
equation provides the energy first integral of motion of the following form: 

u 2 c 2 X 2 

The conserved specific energy £ does not depend on the flow configuration for obvious reason. 
Since r dependence of H varies for different flow geometries, the integral solution of the 
continuity equation, which is another first integral of motion and is referred to as the mass 
accretion rate, will be different for three different accretion configurations. Expressions for 
the mass accretion rate can be obtained as: 

M CH = purH c (4.5a) 
M CM Spur 2 (4.5b) 

IT 3 _I 

M YE = \J-uc s pr2 ($') 2 ( 4>5c ) 

where the subscript CH,CM and VE stands for the flow with constant height (CH), in conical 
model (CM), and in vertical equilibrium (VE), and implies that the respective algebraic 
equations are to be solved to obtain the critical point for the corresponding flow geometries. 
The quantity H c is the constant disc height and 6 is the solid angle sustained by the flow. 
The mass accretion rate for flow in vertical equilibrium not only depends on the matter 
geometry (through the radial dependence of H), but also the information about the space 
time geometry is encrypted in Mve through the explicit appearance of the derivative of the 
pseudo-Schwarzschild potential. One defines the entropy accretion rate as ([14, 36]): 

M M-f^K^ (4.6) 

Substitution of M from eq. (4.5) in the above equation provides the expression for M in 
terms of the adiabatic sound speed, radial distance, and the dynamical flow velocity. The 
space gradient of the sound speed as well as the dynamical velocity for various flow geometries 
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can be obtained by differentiating eq. (4.4) and eq. (4.6): 
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The critical point conditions may be obtained by simultaneously making the 
and the denominator of eq. (4.8) vanish, and the aforementioned critical point 
may thus be expressed as: 
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(4.9b) 
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(4.10a) 
(4.10b) 
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(4.11a) 
(4.11b) 



The critical point conditions for the constant height flow, the conical model flow and 
flow in vertical equilibrium are stated in eq. (4.9), (4.10) and (4.11) respectively. Linearizing 
the Euler and the continuity equation for accretion in hydrostatic equilibrium in the vertical 
direction, one can show that the linear perturbation propagates with the speed ^j^c s 

instead of c s . We thus define J j^c s to be the 'effective' sound speed for such flow. This 
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happens because for flow in vertical equilibrium the expression for the flow thickness contains 
the adiabatic sound speed c s — jp/p. 

For a set of fixed values of [£, A, 7], the location of the critical point for a particular 
flow model can be obtained by substituting the corresponding critical point condition (as 
expressed in eq. (4.9 - 4.11) in the energy first integral eq. (4.4) for a particular pseudo- 
Schwarzschild black hole potential. Once these expressions are substituted, the energy first 
integral becomes an algebraic expression of r c . Exact value of r c for the constant height flow, 
flow in conical model, and in hydrostatic equilibrium can thus be obtained by solving the 
following equations 



The exact location of r c can be evaluated once the astrophysically relevant range of 
[£, A, 7] can be realized. One can argue [13] that the relevant values in the parameter space 
{£, A, 7} can be set as [l < £ < 2, < A<2, 4/3<7<5/3] . 

A solution of eq. (4.12) may exhibit either one (saddle type), or three (one centre type 
flanked by two saddle type) critical points depending on the chosen set of parameters [£, A, 7] 
used. Certain [£,A, 7] mc C [£, A, 7] thus provides the multi critically in accretion solutions, 
where the subscript 'mc' stands for 'multi critical'. The acoustic horizon are thus the col- 
lection of the 'sonic' points where the radial Mach number becomes unity. Such a horizon 
is located on the combined integral solution of eq. (4.7) and eq. (4.8). For inviscid flow, 
a physically acceptable transonic solution which passes through a saddle type sonic point 
can be realized. Such a solution would be an example which confirms the hypothesis that 
every saddle type critical point is accompanied by its sonic point but no centre type critical 
point has its sonic counterpart. For an axisymmetric configuration, in all three geometries 
discussed in this work, a multi-critical flow is thus a theoretical abstraction where three 
critical points (out of which one is always a centre type, through which the integral solu- 
tion can never pass) are obtained as a mathematical solution of the energy conservation 
equation (through the critical point condition), whereas a multi-transonic flow is a realistic 
configuration where accretion solution passes through two different saddle type sonic points. 
One should, however, note that a smooth accretion solution can never encounter more than 
one regular sonic points, hence no continuous transonic solution exists which passes through 
two different acoustic horizons. The only way the multi transonicity could be realized as a 
combination of two different otherwise smooth solutions passing through two different saddle 
type critical (and hence sonic) points and are connected to each other through a discontin- 
uous shock transition. Such a shock has to be stationary and will be located in between 
two sonic points. For a specific [£, A, 7] No Shock C [£, A, 7] mc , three critical points (two saddle 
embracing a centre one) are routinely obtained but no stationary shock forms. Hence no 
multi transonicity is observed even if the flow is multi-critical, and real physical accretion 
solution can have access only to the outer type saddle point out of the two. Thus multi 
critical accretion and multi transonic accretion are not topologically isomorphic in general. 
A true multi-transonic flow can only be realized for [£, A, 7]shock ^ ^> 7] mc , if the criteria 
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for forming a standing shock are met (for details about such shock formation and related 
multi-transonic shocked flow topologies see [17, 27, 37]). For mono-transonic flow, one can 
have only one acoustic horizon on which the related surface gravity may be evaluated. For 
multi-transonic shocked accretion, however, one can have two acoustic horizons (at the inner 
and the outer saddle type critical point) and can calculate the corresponding two different 
values of the acoustic surface gravity. We show this in subsequent sections. 

Once the critical point is located, the critical derivatives of the sound speed (^r) r and of 
the flow velocity (y) , evaluated at the critical point r c (which coincides with the location 
of the acoustic horizon can be obtained for various flow models by applying L' Hospital's 
rule to the numerator and the denominator of eq. (4.8) 
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The quantity u c in (4.10) may be substitutted from eq. (4.11). 

The acoustic surface gravity k as defined in eq. (3.6) may now be evaluated for various 
space time geometries for adiabatic accretion. The location of the acoustic horizon (the 
critical point r c ) and [u,c s ,dc s /dr,du/dr] r can be evaluated as a function of the initial 
boundary conditions as defined by the parameters [£, A, 7] for the adiabatic flow and [T,\] 
for the isothermal flow (the details of the calculation of k for the isothermal accretion will be 
presented in the next section) for a fixed flow geometry in all four pseudo potentials as well 
as under the influence of a particular pseudo potentials in all three different flow geometries 
as considered in this work. 
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4.2 Isothermal Accretion 



For isothermal flow, the integral solution of the time independent Euler equation provides 
the following first integral of motion 



u 



+ c 2 s lnp 



2r 2 



+ $(r) = Constant 



(4.14) 



Obviously, this constant of motion can not be identified with the specific energy of the flow. 
The isothermal sound speed is proportional to T2. The mass accretion rate, another first 
integral of motion of the accreting system of aforementioned kind, may be obtained for three 
different flow geometries as 
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The space gradient of the velocities for these three models comes out to be 
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which provides the following critical point conditions 
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(4.15c) 
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(4.16b) 
(4.16c) 
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for flow with constant thickness (eq. 4.17), in conical equilibrium (eq. 4.18) and for flow in 
hydrostatic equilibrium in the vertical direction (eq. 4.19), respectively. A two parameter 
input [T, A] (T being the isothermal flow temperature), can solve eq. (4.17) - (4.19) to obtain 
the location of the acoustic horizon for three different flow configurations as mentioned above. 
The critical space gradient of the dynamical velocity as evaluated on the acoustic horizon are 
given by 



duY so 
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(4.20b) 



and 
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Hence the acoustic surface gravity for isothermal accretion can be evaluated for three different 
flow models as a function of only two parameters, namely, the flow angular momentum A and 
the isothermal flow temperature T. 

5 Analytical calculation of the acoustic surface gravity and the corre- 
sponding analogue Hawking temperature 

In this section we calculate the surface gravity k for a fluid gravitating in the Paczyhski and 
Wiita (1980) pseudo-Schwarzschild potential <E>i = — 2 ( r -i) f° r now with constant height, in 
conical shape and in hydrostatic equilibrium in the vertical direction, for both the adiabatic 
as well as the isothermal accretion, and will study the variation of k for three different flow 
geometries used. In addition, we calculate k for the same set of initial boundary conditions 
for both the adiabatic and the isothermal accretion under the influence of all four pseudo 
Schwarzschild potentials as defined by (4.2), for a flow in any of the three geometries men- 
tioned before. Studying the acoustic surface gravity k as a function of <E> provides information 
about the dependence of k on the background space time geometry. 

5.1 Adiabatic Accretion 

In subsequent sections, we will calculate k for three different models using Paczyhski and 
Wiita potential [32] for adiabatic accretion. 

5.1.1 Accretion flow with constant thickness 

We start with the simplest flow configuration - axisymmetric flow with constant thickness. 
For such flow, the space gradient of the speed of sound and the flow velocity can be computed 



as: 
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The corresponding critical point conditions can thus be obtained as: 




(5.2) 
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By substituting the above condition into the equation for the energy first integral (4.4), a 
fourth degree polynomial of r c of the following form can be obtained: 



A + r ir 3 c + r 2 r 2 c + r 3 r c + r 4 = o 



(5.3) 
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The location of the acoustic horizon in terms of [5, A, 7] can be obtained analytically by by 
solving the algebraic equation (5.3) for r c using the Ferrari's method (for the details of the 
Ferrari's method and its use in classical algebra, see, e.g., [38]). Then, from eq. (5.2) one 
finds the flow velocity and the sound speed for each solution r c . The critical space gradient of 
the flow velocity and the sound speed evaluated on the acoustic horizon can then be obtained 
by applying l'Hospital's rule on the numerator and the denominator of du/dr in eq. (5.1), 
and then by substituting the value of (du/dr) rc in the expression of dc s /dr in eq. (5.1) on 
the acoustic horizon: 



( d£s 
\ dr 



du 
dr 



r c 



Uc 



7 



1 + 7 

1-7 
1 + 7 



' (1-7) 2 
(1+7) 



u 2 (l - 3 7 ) 



+ 



1 



3A 2 



4(1 + 7 )r 2 4(r c - l) 3 Arf 



(5.4) 



l u 2 (l - 3 7 ) 1 3A 2 

r 2 (l + 7 ) 2 + (l + 7)(r c -l) 3 r|(l + 7 ) 



Note that the quantities r c , c Sl dc s /dr and du/dr evaluated at the acoustic horizon are 
expressed in terms of elementary functions of £, A, and 7. Hence, the surface gravity k can 
be calculated analytically as a function of [£, A, 7] since 



^CH — CcH ( r i C Si ~T~ ' 



dc s du \ 
dr ' dr J 



(5.5) 



as is obvious from eq. (3.6). 



5.1.2 Conical Model 

For conical flow, the space gradient of c s and u can be obtained as 



dc s 
dr 



c s (l- 7 ) 



du_ u _ 
dr 

Hence the critical point conditions becomes 

(c a ) = (u) = >j 



— + K 



1 du 2 

u dr r 
1 



(5.6) 



2(r-l) 2 



(u 2 - C 2 ) 



4(r c 



A 2 



l) 2 2r 2 



(5.7) 
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The corresponding fourth degree polynomial in r c can be expressed as 

rj + rv c 3 + r 2 r c 2 + r 3 r c + r 4 = o 

where 

(3 7 -5)-16g( 7 -l) p 2( 7 -l)(2£:-l)-A 2 (7-3) 

1 1 = T7T7 77 , 1 2 = 



(5. 



85( 7 - 1) 



45( 7 - 1) 

A 2 ( 7 -3) A 2 (3- 7 ) 
:> 1 4 ~ 



25( 7 -l)' 45( 7 -l) 
The critical gradient of the sound speed and the flow velocity can be obtained as 

2u r . 



( dcs 
I dr 



Ml-7) 4 A2(2-r c )(l- 7 ) 2 (l-7) 2 [r c (3-2 7 ) + (27-lj] 
VI + 7/ V r c(l + 7) 2 2r c 4 (l + 7 ) 8r c (l + 7 )(r c -l) 3 1 
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dr 



2uc A~7 
r c \l + 7 



4ix§ A - 7 V 2A 2 (2 - r c ) r c (3 - 2 7 ) + (2 7 - 1 



1 + 7 



+ 



(l + 7 )r c 4 2r c (l + 7 )(r c -l) 



Using eq. (5.7) - (5.10) the acoustic surface gravity for the conical flow 



dc s du 
k cm — Ccm I r i c s i ~r~ 5 ~r~ 



(5.10) 



can thus be calculated analytically as a function of [£, A, 7]. 



5.1.3 Flow in hudrostatic equilibrium in vertical direction 

For flow in hydrostatic equilibrium in the vertical direction, the velocity gradients and the 
corresponding critical point conditions become 



dc s 
dr 



7-1 
7 + 1 



-1 du 



5r-3 



u dr 2r(r — 1) 



(5.11) 
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The corresponding fourth degree polynomial in r c can be expressed as 

r 4 c + T ir 3 c + T 2 r 2 + T 3 r c + T 4 = 



Ti = 



c 1 J-' c 1 4' c 

where 

5 - 165 - 65 - 3 + ?4A 2 



(5.12) 



(5.13) 
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r, r 4 
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° 10(7-1)5' 10(7-1)5 
The critical gradient of the flow velocity can be found as 



du\ 



-p - a//3 2 - 4aS 
2a 



(5.14) 
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where 



A 2 ( 7 + l)(r c -2) 
ri(r c - 1) 



7+1 



A 2 (5r c -3)( 7 -l) 



2r c (r c - l) 3 
5( 7 + 1) 



a = 47 
2( 7 -l)(5r c -3K 
r c {r c ~ 1) 
, (5r c -3)( 7 -l) 
2r c (r c - l) 3 
5A 2 ( 7 + 1) n 



+ 



r*(r c -l) 2(r c -l) 2 (5r c -3) r c 3 (5r c - 3) 

Hence the critical gradient of the speed of sound can be found as 



(5.15) 
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'(7 + l)(r c 3 
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-1) 2 A 2 ) 


2r c 2 (r c - 


-l)(5r c 


-3) 



7-1 
7 + 1 



-1 fdu\ 
u c \drj, 



5r c - 3 
2r c (r c - 1) 



(5.16) 



where (^) r is to be substituted from eq. (5.14). 
5.2 Isothermal accretion 

For isothermal flow under the influence of the Paczyhski and Wiita (1980) potential, specific 
energy does not remain one of the first integrals of motion any more. The mass accretion 
rate, however, still remains a constant of motion. Since the temperature is constant, the 
value of the isothermal sound speed c s = \/ i^b / {i^^h)T^ is position independent and hence 
dc s /dr — identically. Mach number profile for the isothermal accretion is thus found to be 
a scaled down version of the dynamical velocity profile. The stationary solution is completely 
characterized by two parameters [T, A] , T being the isothermal flow temperature. 

5.2.1 Accretion flow with constant thickness 

For constant thickness flow, the velocity gradient 



du 
dr 



u 



+ 4 



2(r-l) 2 



{u 2 - eg) 



(5.17) 



provides the critical point condition as 

0)r c = (°s)r c = V^B/(^rn H )T^ 
the corresponding fourth degree polynomial is 



A 2 



L2(r c - l) 2 



rt + rv c 3 + r 2 r 2 + r 3 r c + r 4 = o 



where 
-2A 2 _ A 2 



(5.18) 
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1 A 2 

ri = - 2 "2c!' r2 = 1+ c!' r3 = 

The critical flow velocity gradient becomes 
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(dr 



r c + l 



A 2 



4r c (r c 



(5.20) 
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Although both c s and du/dr evaluated at the acoustic horizon depend only on the angular 
momentum of the flow and not on the flow temperature, the location of the acoustic horizon 
itself (the critical point r c ) is a function of both T and A, hence 



r c — 

SCH \ 1 1 ^ 



ASO /-ISO 



(5.21) 



can be calculated analytically (for all three different flow models considered here, as we will 
see in subsequent sections) as a function of only two accretion parameters [T, A] . 

5.2.2 Conical Model 

For conical flow the velocity gradient 



du_ v _ 
dr 



2c? , A 2 



2(r-l) 2 



(u 2 - eg) 



(5.22) 



provides the critical point condition as 

(«)r c = ( c »)r c = V K B/(^m H )T^ = 
The corresponding polynomial in r c becomes 



A 2 



4(r c - l) 2 



r 4 c + IVl? + T 2 r 2 c + F 3 r c + T 4 = 



(5.23) 



(5.24) 



" 2 "4c|' r2 - 1+ 2^' r3 



where 
c 2 ' U 2c 2 



and the critical gradient of the dynamical velocity is thus 



/ du 
(d^ 



r c + l 



A 2 



4r c (r c - l) 3 r c 4 



(5.25) 



which is identical to that obtained for a flow with constant thickness, see eq. (5.20). The 
corresponding acoustic surface gravity 



„lSO A1SO 

V CM SCM 



s ' dr J 



(5.26) 



can thus be calculated analytically as a function of [T, A]. 

5.2.3 Accretion in hydrostatic equilibrium in vertical direction 

For flow in hydrostatic equilibrium in the vertical direction, the corresponding quantities are 

" <2(5r-3) + 



du 
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2r(r-l) 



2(r- 1) 2 



(u 2 - C 2 ) 
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' 2r e (r c -l) 
(5r c - 3) 



1 A 2 
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r 4 c + 1^ + T 2 r c 2 + T 3 r c + T 4 = 

where 
-4A 2 r4 _ 2A^ 
5c? ' 5c? 



1 + 8 C 2 ^ 3c 2 + 2A 2 

2 — 1 1 3 
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du\ 
to). 



2r c (r c - l)(5r c - 3) 



5r 2 -3 2A 2 (5r 2 - 9r c + 3) 



2(r c - l) 2 



and the corresponding acoustic surface gravity 



^ve = We \ T i c si J 



^lSO A1SO 



(5.30) 



(5.31) 



can be evaluated accordingly. 



6 Dependence of the acoustic surface gravity on the flow geometry and 
initial boundary conditions 



Figure 1 shows the variation of the acoustic surface gravity z^asa function of the specific 
angular momentum A for monotransonic polytropic accretion. The specific energy £ and 
the polytropic index 7 have been kept constant at values [£ — 0.06,7 — 1-333]. The range 
of A for which k has been calculated for a particular flow model constructs a subset in the 
[£, A, 7] parameter space for which the representative flow models produces mono-transonic 
accretion for a fixed set of values of [£]. Alternative ranges for A for other similar subsets 
of [£,A, 7] parameter space may also be considered to study the 'ft — A' profile for other 




Figure 1. Variation of the acoustic surface gravity with the specific angular momentum of the flow 
for monotransonic polytropic accretion characterized by £ = 0.06, 7 = 1.333 for three different flow 
geometries - namely, for flow in hydrostatic equilibrium in the vertical direction (solid green line), 
constant thickness flow (long dashed black line), and for the conical model (dotted red line). 
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Figure 2. Variation of the acoustic surface gravity with the specific energy of the flow for monotran- 
sonic poly tropic accretion characterized by A = 1.835, 7 = 1.333 for three different flow geometries - 
namely, for flow in hydrostatic equilibrium in the vertical direction (solid green line), constant thick- 
ness flow (long dashed black line), and for the conical model (dotted red line). Only the common 
range of specific energy for which all the flow configurations produce adiabatic monotransonic flow 
has been considered in the figure, see text for further detail. 




Figure 3. Variation of the acoustic surface gravity with the specific energy of the flow for monotran- 
sonic poly tropic accretion characterized by A = 1.85, £ = 0.06 for three different flow geometries 
- namely, for flow in hydrostatic equilibrium in the vertical direction (solid green line), constant 
thickness flow (long dashed black line), and for the conical model (dotted red line). 



fixed values of [£,7]. The solid line (green coloured in the online version) represents the L n 
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Figure 4. Variation of the ratio of the acoustic surface gravity evaluated at the inner acoustic horizon 
to the acoustic surface gravity evaluated at the outer acoustic horizon with the range of the specific 
angular momentum (for a fixed value of [£ = 0.0002,7 = 4/3]) of the multi-transonic flow for which 
a stationary shock may form. 



- A' variation profile for mono-transonic accretion in vertical equilibrium (VE), whereas the 
dotted (red coloured in the online version) and the long dashed (black coloured in the online 
version) curves represents such dependence for flow with constant height (CH) as well as for 
the conical flow (CM) respectively. 

It is usually observed that k varies with A non linearly and non monotonically. For rela- 
tively lower values of the specific angular momentum, the acoustic surface gravity correlates 
with the specific angular momentum and attains a peak characterized by an unique value 
of the specific angular momentum denoted by A max , and subsequently falls of with A for 
A > A max . A max is different for different flow models and one observes that 

\CM \ VE \ CH / fi n 
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Figure 5. Variation of the acoustic surface gravity with the specific angular momentum of the flow 
for monotransonic isothermal accretion characterized by the isothermal flow temperature Xio = 22 for 
three different flow geometries - namely, for flow in hydrostatic equilibrium in the vertical direction 
(solid green line), constant thickness flow (long dashed black line), and for the conical model (dotted 
red line). 

For a set of fixed values of [£, A] , A max for a particular flow model can be calculated completely 
analytically. We illustrate such procedure for constant thickness flow. Similar procedure may 



0.12 




Figure 6. Variation of the acoustic surface gravity with the isothermal flow temperature of the flow 
for monotransonic isothermal accretion characterized by A = 1.8 for three different flow geometries 
- namely, for flow in hydrostatic equilibrium in the vertical direction (solid green line), constant 
thickness flow (long dashed black line), and for the conical model (dotted red line). 
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be followed to calculate A max for two other flow geometries. 

For constant thickness flow, eq. (5.5) provides the dependence of ft on the critical points, 
the acoustic speed and its space gradient, and on the space gradient of the flow velocity 
itself, everything evaluated at the acoustic horizon (the critical point). (c s ) , (dc s /dr) rc and 
(du/dr) r can be expressed as a function of r c and [£,A,7] using eq. (5.2) and eq. (5.5), 
respectively. The critical point r c itself can be computed in terms of [£, A, 7] by solving the 
polynomial in r c as represented through eq. (5.3). Hence ft for the constant thickness flow, 
can be specified in terms of [£, A, 7]. Analytical expression for ftcH=^CH [£•> A, 7] 

can thus be maximized with respect to the specific angular momentum and the corre- 
sponding A max can thus be obtained. 

From figure 1, one should note that for the common range of the specific angular mo- 
mentum for which all three flow models will produce mono transonic flow for a fixed set of 
value of [£,A], does not allow to explore the complete non monotonic 'ft - A' profile for all 
the three flow models simultaneously. For example, for the common range of A as shown in 
the figure for which all three flow geometries produces the monotransonic accretion, ft for 
constant height flow as well as for conical model accretion will anti correlates with A, whereas 
for full range of allowed A to form mono transonic accretion at individual level, 'ft - A' profile 
for both the aforementioned flows exhibits a maximum. 

Similar features are observed for the 'ft - £ ' profile, which has been shown in the figure 2, 
where the ft vs £ variation (for a fixed set of [A, 7]) for constant height flow apparently 
shows that ftcH correlates with £ , whereas ^cm and ftyE anti-correlates with £. Such trend, 
however, does not provide the complete information about the 'ft - £' profile in general since 
the range of £ for which the figure is drawn is actually taken from the common region of 
[£,A, 7] space for which all three flow geometries produces mono transonic accretion for a 
fixed value of [A, 7]. If one allows ft to vary with £ for the entire range of the specific energy 
for which a particular flow model provides the mono transonic accretion for a fixed value 
of [A, 7], ft vs £ profile would have a non monotonic behaviour with a corresponding £ max 
separately for every flow configuration. The corresponding £ max for every flow model could 
then be estimated by maximizing the expression for the acoustic surface gravity with respect 
to £ by keeping [A, 7] constant. One thus understands that the 'anti correlating' ft vs £ 
profile for the conical flow (represented by the dotted red curve) as well as for flow in the 
hydrostatic equilibrium along the vertical direction (solid green curve), respectively, are thus 
the post-peak {£ > £ max ) descending part of the complete non monotonic 'ft - £' profiles 
for the corresponding flow configuration. Similarly, the 'correlating' ft - £ profile for the 
constant height flow (represented by the long dashed black curve) is actually the pre-peak 
{£ < £ m ax) ascending part of the complete non monotonic ft - £ profile for the corresponding 
flow geometry. Hence f^ax ls largest among all the values of £ max corresponding to all three 
different flow configurations. 

Figure 3 shows the ft - 7 profile (for a fixed set of [£ , A]). It is obvious from the figure that 
7max nas the largest value among all three values of 7 max corresponding to three different flow 
configurations. £ max and 7 max for any particular model can also be estimated by maximizing 
ft with respect to the respective parameters. 

One thus understands that for monotransonic adiabatic accretion characterized by a fixed 
set of values of [£,A, 7], the analogue surface gravity for three different flow configurations 
exhibit the following trend 

adia ^ adia ^ adia (n r>\ 

ft CM ^ ^VE ^ ^CH 
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Hence for the adiabatic mono transonic accretion, flow in conical shape produces the largest 
value of the analogue Hawking like temperature whereas the flow with constant thickness 
provides the lowest value of Tah f° r the same set of initial boundary conditions. 

As has been mentioned in section 4.1, multi-transonic accretion with stationary shock 
may be realized for all three different flow configurations for adiabatic as well as for isother- 
mal accretion, studied here in this work. For such flow topologies, two black hole type 
acoustic horizons form at the inner and the outer saddle type sonic points. The correspond- 
ing acoustic surface gravity K\ n and K out can be evaluated at the inner and the outer sonic 
points respectively. We define 

Kin (n o\ 

K io (6.3) 

^out 

It has been observed that the overall — [£, A, 7]) as well as the (n ut — [£, A, 7]) profiles 
are qualitatively similar with the (k — [£, A, 7]) profile for all three flow configurations, where 
k is the acoustic surface gravity evaluated for the mono transonic accretion. For all three flow 
geometries considered in this work, the value of K out is, however, several orders of magnitude 
less than that of the corresponding K{ n for same set of initial boundary conditions. This is 
because the outer acoustic horizons form at a large distance away from the black hole event 
horizon compared to the location of the inner acoustic horizon with respect to the black hole 
event horizon. For a typical set of the initial boundary conditions, the inner acoustic horizon 
may form 1.5-5 Schwarzschild radius away from the black hole event horizon whereas the 
outer acoustic horizon may be located at 10 3 - 10 6 Schwarzschild radius away, or even more, 
from the black hole event horizon. Weak gravity at such a large distance (where the outer 
sonic horizons form) restrict the acoustic surface gravity to possess such a small numerical 
value. Same argument applies for the analogue Hawking temperatures evaluated at the inner 
and the outer acoustic horizons as well. 

In figure 4, we plot fcj%o as a function of A for the constant thickness flow (uppermost 
panel), flow in hydrostatic equilibrium along the vertical direction (mid panel) and for the 
conical flow (lowermost panel). The range of A used in this figure corresponds to the common 
value of the specific angular momentum for which shock forms for multi-transonic accretion 
in all three flow geometries for a fixed value of [£ , A] . The common range of A is chosen in such 
a way so that the inner acoustic horizon forms at a distance larger than two Schwarzschild 
radius from the black hole event horizon. Since we use pseudo Schwarzschild potentials which 
are relatively less reliable in simulating the general relativistic space time extremely close to 
the black hole event horizon, we prefer to confine our attention to the mono transonic as well 
as the multi transonic flow topologies for which the acoustic horizon does not form at a very 
close proximity to the black hole event horizon. From the figure it is evident that 

\CM . \VE . \ CH / fi 4 x 

A max ^ A max ^ A max V u / 

where A max in eq. (6.4) is the value of A for which the non monotonic K{ - A profile attains its 
maximum. Interestingly enough, eq. (6.1) is identical with eq. (6.4), indicating the fact that 
the dependence of the acoustic surface gravity on initial boundary conditions is similar for 
both the mono as well as for the multi-transonic shocked accretion in all three flow geometries 
considered here in this work. Once again, the K{ - A dependence does not exhibit the peaked 
non monotonic profile since the common range corresponding to the shock forming A is 
not sufficient to provide the required span for the K{ - A variation for any particular flow 
geometry for the entire range of A for which the corresponding flow model produces shocked 
multi-transonic accretion for a fixed value of [£, A]. 
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Figure 5 demonstrates the dependence of the acoustic surface gravity on specific angular 
momentum for mono transonic isothermal accretion for three different flow models. It is 
observed that 

A CM \ ve \ CH / fi ,\ 

In figure 6, we plot the dependence of k on the isothermal flow temperature T (scaled by a 
factor of 10 10 degree Kelvin - T w =Tx 10" 10 ). We obtain 

r CM r VE r CH , g 6) 

-'-max ^ max ^ max V w " w / 

The value of A max and T max for the monotransonic isothermal flow may be estimated for all 
three flow geometries in a way similar to what has been accomplished for the polytropic flow. 



7 Discussion 

Axisymmetric accretion onto a non-rotating astrophysical black hole under the influence of 
pseudo-Schwarzschild potentials are natural example of classical analogue systems found in 
the universe. The corresponding acoustic geometry may be studied for three different flow 
configurations, viz, accretion with constant flow thickness, the conical flow, and accretion 
disc in hydrostatic equilibrium in the vertical direction. For each background flow geometry, 
all together eight different configurations for the acoustic geometry - adiabatic as well as 
isothermal accretion for four different pseudo-potentials, may be studied. For any specific 
pseudo potential, six different configurations of the acoustic flow geometry - adiabatic as well 
as isothermal accretion in three different flow geometries, may be studied. For any geometric 
configuration of the adiabatic flow under the influence of a particular potential, the three 
initial parameters, viz., the specific energy S , the specific angular momentum A, and the 
adiabatic index of the flow 7 can completely specify the corresponding acoustic geometry. 
Similarly for the isothermal flow in any potential and with any flow geometry, the two initial 
parameters, viz., the constant flow temperature T and the specific angular momentum A, can 
completely specify the corresponding acoustic geometry. 

Among six possible flow configurations for any particular pseudo potential used, only the 
adiabatic flow in vertical equilibrium exhibits an 'effective' sound speed which is a scaled ver- 
sion of the adiabatic sound speed with a 7 dependent scaling constant. The scaling constant 
becomes unity for isothermal flow. The reason is that for the accretion in vertical equilib- 
rium the flow thickness is a function of the adiabatic sound speed, as well as a function of the 
space derivative of the pseudo-potential used, since the expression for such flow thickness is 
obtained by balancing the pressure gradient with the relevant component of the gravitational 
force. One should, however, bear in mind that the corresponding expressions for the flow 
thickness in all three flow geometries used are derived using a set of idealized assumptions. In 
principle, a more realistic derivation of the flow thickness may be worked out by employing 
the non-LTE radiative transfer (see [39, 40] ) or by taking recourse to the Grad-Shafranov 
equations for MHD flow (see [41-43] ). 

For multi-transonic flow, two acoustic black holes are formed at two regular saddle type 
sonic points, whereas the acoustic white hope forms at the shock location, in agreement with 
the results obtained by [19] and [8] . The acoustic surface gravity is formally infinite for the 
acoustic white hole since the flow velocity as well as the sound speed changes discontinuously 
at the shock location, in agreement with the findings of [44] . 

The surface gravity k (or Tah) profile obtained for the multi-transonic flow at the inner 
acoustic horizon is similar to the k profile for mono-transonic flow. This indicates that 
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irrespective of the phase topology, the surface gravity is basically determined by the physical 
proximity of the acoustic horizon to the black hole event horizon. This is further supported by 
the fact that irrespective of the flow geometry, pseudo potentials or the equation of state used 
to describe the accretion flow, the value of k at the outer acoustic horizon (for multi-transonic 
flow) is much less than that evaluated at the inner acoustic horizon. 

For a fixed set of [£,A, 7] describing the adiabatic accretion as well as for a fixed set 
of [T, A] describing the isothermal accretion, the conical flow produces the largest whereas 
the flow with constant thickness the smallest surface gravity and the analogue Hawking 
temperature. The accretion flow in the hydrostatic equilibrium in the vertical direction 
provides the value of n and Tah in between the respective values of k and Tah for the 
conical and the constant thickness flow, respectively. The influence of the flow geometry in 
determining the acoustic surface gravity has thus been successfully investigated in this work. 

One of our major achievements is the computation of the acoustic surface gravity and the 
investigation of its dependence on various flow geometries and on various accretion parameters 
completely analytically.. However, one should note that this has been possible owing to a 
specific feature of the chosen pseudo potential. With the Paczyhski & Wiita (1980) potential 
<E>i = — 2 (r-i) ? the energy first integral for the adiabatic accretion (see eq. (4.12) as well as 
the critical point condition for the isothermal flow (see eq. (5.5 - 5.7)) can be recast into a 
fourth degree polynomial in r c (equivalent ly, in 77J in an exactly solvable form. This has not 
been possible with other pseudo potentials, for which such an exactly solvable polynomial in 
rh can not be constructed. With other potentials, however, the surface gravity k can still 
be studied as a function of various flow parameters and flow geometries, with the help of 
numerical methods. Remarkably, it has been demonstrated in the literature ([16, 17] and 
references therein) that out of the four pseudo Schwarzsculd potentials as shown in eq. (5.9), 
the potential $1 mimics the Schwarzschild space time most efficiently in constructing the 
integral flow solution for transonic accretion. 

However, even if one can not employ a complete analytical calculation and even if $1 is 
the most suitable potential to mimic a Schwarzschild space time, our generalized formalism 
for the evaluation of [u,c s ,dc s /dr,du/dr] rh and the corresponding values of k and Tah in 
terms of a general $ is still important in the following sense. There exists a possibility that a 
new form of <E> more effective than $1 will be suggested in the future in order to approximate 
the Schwarzschild space time in constructing the integral accretion solutions for the transonic 
flow. In such a case, if a general model is capable of computing the acoustic surface gravity 
(and hence the analogue Hawking temperature) in a way presented in this work, then this 
model will be able to readily accommodate that novel form of the pseudo-Schwarzschild 
potential with no need to significantly change the fundamental structure of the formulation 
and the solution scheme. In this case one need not worry about providing any new unique 
scheme valid exclusively only for a particular form of pseudo potential. 

The methodology developed in this paper can also be used to construct the relevant 
acoustic geometry for the equatorial slice of the accretion flow under the influence of various 
pseudo-Kerr potentials as well to study the variation of the surface gravity with the black 
hole spin. Such work is in progress and will be reported elsewhere. 

However, one should be cautious when using the pseudo potentials because none of the 
potentials discussed here can be directly derived from the Einstein equations. These poten- 
tials are used to obtain more accurate correction terms over and above the pure Newtonian 
results. Hence any 'radically new' results obtained using these potentials should be cross 
checked very carefully against general relativity. Besides, one should bear in mind that these 
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potentials are not too reliable for modeling the space time in the very close neighborhood of 
the event horizon since strong gravity effects dominate in this region. Hence our formalism 
may not be quite realistic if the acoustic horizon forms very close to the black hole event 
horizon. We thus consider only those initial boundary conditions for which r\ l > 2 ensuring 
our findings to be trustworthy. 
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